Memory Elements: A Paradigm Shift in 
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Abstract: Meminductors and memcapacitors do not allow a Lagrangian formulation in the 
classical sense since these elements are nonconservative in nature and the associated energies 
are not state functions. To circumvent this problem, a different configuration space is considered 
that, instead of the usual loop charges, consist of time-integrated loop charges. As a result, 
the corresponding Euler-Lagrange equations provide a set of integrated Kirchhoff voltage laws 
in terms of the element fluxes. Memristive losses can be included via a second scalar function 
that has the dimension of action. A dual variational principle follows by considering variations 
of the integrated node fluxes and yields a set of integrated Kirchhoff current laws in terms 
of the element charges. Although time-integrated charge is a somewhat unusual quantity in 
circuit theory, it may be considered as the electrical analogue of a mechanical quantity called 
absement. Based on this analogy, simple mechanical devices are presented that can serve as 
didactic examples to explain memristive, meminductive, and memcapacitive behavior. 
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1. INTRODUCTION 

Memristors, meminductors, and memcapacitors constitu- 
tive an increasingly important class of two-terminal circuit 
elements whose resistance, inductance, and capacitance re- 
tain memory of the past states through which the elements 
have evolved. All three elements are nonlinear and can be 
identified by their pinched hysteresis loop in the voltage 
versus current plane, current versus flux plane, and voltage 
versus charge plane, respectively. While there are many 
discovered experimental realizations and applications of 
systems that exhibit memristive behavior, ranging from 
applications in non-volatile nano memory to intelligent 
machines with learning and adaptive capabilities, the num- 
ber of systems showing memcapacitive and meminductive 
behavior is still somewhat limited. Nevertheless, several 
applications of these concepts are foreseen in the field 
of logic and arithmetic operations using memristive and 
memcapacitive devices, and field-programmable quantum 
computation using meminductive and memcapacitive de- 
vices. See Di Ventra et al. (2009) for an introduction into 
the memory elements, and Pershin and Di Ventra (2011) 
for a fairly complete overview of the current state-of-the- 
art and an extensive list of references. 

In the past century, a significant amount of research has 
been devoted to the description and analysis of conven- 
tional electrical circuits using Lagrangian and Hamilto- 
nian methods; see Jeltsema and Scherpen (2009) and the 
references therein. Since these classical frameworks are 
important not just for their broad range of applications, 
but also for their role in advancing deep understanding of 
physics, the next step is to consider circuits made from 
memristors, meminductors, and memcapacitors. 



The remainder of the paper is organized as follows. In 
Section 2, the mathematical properties of the memris- 
tor, meminductor, and memcapacitor are discussed. The 
structure of the Lagrangian equations for a large class of 
conventional circuits is briefly reviewed in Section 3 and it 
will be shown that meminductors and memcapacitors do 
not fit into this framework a priori since these elements 
are nonconservative in nature and the associated energies 
are not state functions. To circumvent this problem, in 
Section 4 a different configuration space is considered that, 
instead of the usual loop charges, consist of time- integrated 
loop charges. The Lagrangian is defined by the difference 
between two novel state functions in a fashion similar 
to the usual magnetic co-energy minus electric energy 
setup, but having the dimensions of energy times time- 
squared which, in turn, is equivalent to action times time. 
As a result, the corresponding Euler-Lagrange equations 
provide a set of integrated Kirchhoff voltage laws (iKVL's) 
in terms of the element fluxes. Memristive and resistive 
losses can be included via the introduction of a second 
scalar function that has the dimension of action. Further- 
more, a dual variational principle follows by considering 
variations of the integrated node fluxes and yields a set 
of integrated Kirchhoff current laws (iKCL's) in terms 
of the element charges. Although time-integrated charge, 
which in SI units is measured in Coulomb times seconds, 
is a somewhat unusual quantity in circuit theory, it may 
be considered as the electrical analogue of a mechanical 
quantity called absement. Based on this analogy, in Sec- 
tion 5, simple mechanical devices are presented that can 
serve as didactic examples to explain meminductive and 
memcapacitive behavior. Finally, the paper is concluded 
with some final remarks in Section 6. 



2. DEFINITION OF MEMORY ELEMENTS 

From a mathematical perspective, the behavior of a two- 
terminal resistor, inductor, and capacitor, whether linear 
or nonlinear, is described by a relationship between two 
of the four basic electrical variables, namely, voltage V, 
current /, charge q, and flux ip, where 

q(t) = f /(r)dr, (1) 



cp(t) = / y(r)dr. 



(2) 



A resistor is described by a constitutive relationship be- 
tween current and voltage; a capacitor by that of voltage 
and charge; and an inductor by that of current and flux 
linkage. 

Based on logical and symmetry reasonings, Chua (1971) 
postulated the existence of a fourth element that is char- 
acterized by a constitutive relationship between the re- 
maining two variables, namely charge and flux. This el- 
ement was coined memristor (a contraction of memory 
and resistor) referring to a resistor with memory. The 
memory aspect stems from the fact that a memristor 're- 
members' the amount of current that has passed through it 
together with the total applied voltage. More specifically, 
if q denotes the charge and ip denotes the flux, then a 
charge-modulated memristor is defined by the constitutive 
relationship ip — <p(q). Since flux is defined by the time 
integral of voltage V (like in Faraday's law), and charge is 
the time integral of current /, or equivalently, V = ip and 
I = q, we obtain 

V = R M (q)I, (3) 
where Rmiq) '■= dip/dq is the incremental memristance. 

Note that (3) is the definition of a memristor in impedance 
form. The admittance form / = Gm(^)^, with incremen- 
tal memductance Gm^) '■= dq/dtp, is obtained by starting 
from a constitutive relationship q = q(ip). 

In addition to the memristor, it is shown in Di Ventra 
et al. (2009) that the memory-effect can be associated to 
inductors and capacitors as well. For that, let a and p 
denote the time-integrals of charge and flux, 

a(t) = [ <z(r)dT, (4) 



P(t) 



</?(r)d7 



(5) 



respectively. Then, a memory inductor, or meminductor 
for short, is a two-terminal element defined by a constitu- 
tive relation p = p(q). Indeed, differentiation of the latter 
with respect to time yields 

ip = L M {q)I- (6) 
Since LM{q) '■= dp/dq relates flux with current, its values 
clearly correspond to the units of inductance. Bearing 
in mind that charge is the time integral of current, the 
memory aspect of a charge-modulated meminductor stems 
from the fact that it 'remembers' the amount of current 
that has passed through it. 

Dually, a flux-modulated memcapacitor is defined by a 
constitutive relationship a — <r(ip), which, after differenti- 



ation with respect to time, yields 

q = C M (<p)V, (7) 

where Cm{<-p) := da/dip denotes the incremental capaci- 
tance. The memory aspect of a flux-modulated memcapac- 
itor stems from the fact that it 'remembers' the amount 
of voltage that has been applied to it. 

A meminductor (resp. memcapacitor) that depends on 
the history of its flux (resp. charge) can be formulated 
by starting from a constitutive relationship of the form 
q = q(p) (resp. ip = ip(a)). 

In the special case that the constitutive relationship of 
a memristor is linear, a memristor becomes an ordinary 
linear resistor. Indeed, in such case (3) reduces to ip = Rq, 
with constant memristance M (the slope of the line), or 
equivalently, V — RI, which precisely equals Ohms law. 
Hence it is not possible to distinguish a two-terminal linear 
memristor from a two-terminal linear resistor. The same 
holds for a linear meminductor and linear memcapacitor, 
where (6) reduces to p = Lq, or equivalently, p> = LI, and 
(7) to d = dp, or equivalently, q — CV, respectively. 

3. SELF-ADJOINTNESS OF CIRCUIT DYNAMICS 

The dynamical behavior of any electrical circuit consisting 
of conventional, possibly nonlinear, resistors, inductors, 
and capacitors is basically determined by three types 
of equations: those arising from Kirchhoff's voltage law 
(KVL), those arising from Kirchhoff's current law (KCL), 
and the constitutive relationships of the elements. In many 
cases this leads to differential equations of the form 1 

Aij(x)x j + Bi{x,x) = 0, (8) 

where i,j = l,...,n, and x € K™ represents a column 
vector of loop charges and node fluxes. 

The system of differential equations (8) allow a Lagrangian 
description if we can find a Lagrangian C(x, x) satisfying 

d dC dC ,-, , . % ... 

-g--^-^A J (x)x> + B i (x,x). (9) 

In mechanics it is known that the existence of a Lagrangian 
C(x,x) relies on the fact that the system of differential 
equations (8) is self-adjoint, which for the present form is 
tantamount to the following set of integrability conditions 
(Santilli, 1978): 

(10) 
(11) 
; , (12) 
(13) 
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3.1 Conventional Conservative Circuits 

From a Lagrangian perspective, it is well-known that 
for a large class of circuits that are made from current- 
controlled inductors and charge-controlled capacitors, and 

1 Throughout the document we adopt the summation convention of 
repeated indices. 



consisting of n loops, the differential equations (8) can be 
written as 

dt dtf dq* ' 1 j 

where q % represents the loop charge associated to the i- 
th loop, and q 1 represents the loop current circulating 
in the z-th loop. For this particular case, the Lagrangian 
C{q 1 q) equals the total magnetic co-energy stored in the 
inductors, T*(q), minus the total electric energy stored in 
the capacitors, U(q). 

On the other hand, if the inductors are flux-controlled and 
the capacitors are voltage-controlled, we need to start from 
a node analysis yielding a so-called co-Lagrange equation 
associated to each node in the network. Hence, if ip> 
represents the potential of the j-th node together with 
its time-integral, the node flux (pP , one obtains 
ddC* dC* 

where the co-Lagrangian C*(<p, ip) equals the total electric 
co-energy stored in the capacitors, U*(ip), minus the total 
magnetic energy stored in the inductors, T(tp). 

As an illustration, consider a circuit consisting of a nonlin- 
ear current-controlled inductor, with constitutive relation 
if = <p(I), and a linear capacitor C as shown in Fig. 1. 

/ V 



:c 



Fig. 1. Circuit with a nonlinear conventional inductor. 

Let q denote the loop charge and q = I the associated loop 
current, then the Lagrangian for the circuit reads 

which, upon substitution into (14), yields the equation of 
motion <p' '(q)q + q/C = 0. Note that the latter constitutes 
the KVL for the circuit. 

Dually, if the nonlinear current-controlled inductor is re- 
placed by a nonlinear flux-controlled inductor, we have to 
consider 

C*(^p)^\cp 2 - J l{<p)l<p, 

resulting in the equation of motion Cip + I (if) = 0, which 
constitutes the KCL for the circuit. 



3.2 Conventional N on- Conservative Circuits 

Dissipation due to conventional resistors in the circuit 
leads to non-conservative dynamics that is not self-adjoint. 
Resistors can therefore not be included using a (standard) 
Lagrangian function. 2 

2 In some cases the dynamics can be rendered self-adjoint by looking 
for an integrating factor (Santilli, 1978). However, the form of the 
differential equations is altered and the Lagrangian does not have 




To circumvent this problem, a so-called content function, 
which is a nonlinear multi-domain generalization of the 
Rayleigh dissipation function, is usually introduced. 

As an example, suppose the inductor in Fig. 1 possesses a 
constant series resistance R. In that case, the equation of 
motion extends to <p'{q)q + Rq + q/C = 0, so that A(q) = 
(f'(q) and B(q,q) = Rq + q/C. For R — 0, the conditions 
(10)-(13) are clearly satisfied, but if R ^ 0, condition (13) 
is violated and no Lagrangian can be associated to the 
circuit. On the other hand, the introduction of a function 
of the form V(q) = \Rq 2 , satisfying 

d_dC dC dV_ {) (w) 
dt dq dq dq 

solves the problem. However, the variational character of 
the dynamics that makes the Lagrangian formalism so 
appealing is clearly lost. 

Note that the existence of a content function relies on the 
fact that the resistors in the circuit are current-controlled. 
The inclusion of voltage-controlled resistors requires the 
introduction of a co-content function T>* (tp) . 



4. CIRCUITS WITH MEMORY ELEMENTS 
4-1 Problem 1: Path-Dependence 

Now consider a circuit consisting of a meminductor (6) 
and a linear capacitor shown in Fig. 2. Application of KVL 
yields ip + q/C — 0, or equivalently, 

^(L M (q)l) = L M (q)i + L' M (q)I 2 + ± = 0. (17) 

In order to derive the dynamics using a Lagrangian for- 
mulation, one is tempted to start from a Lagrangian that 
equals the magnetic co-energy stored in the meminductor 
minus the electric energy stored in the capacitor, i.e., 



£(?»«) = \ L M{q)q 2 - ^q 2 - 



(18) 



However, the latter is clearly not a proper state function 
since it depends on the path q. 




:c 



Fig. 2. Circuit with a meminductor. 

Furthermore, substitution of the Lagrangian (18) into (14) 
yields the equation LM(q)q+ \L' M (q)q 2 + q/C = 0, where, 
in comparison to (17), we observe the appearance of an 
erroneous factor |. 



the usual interpretation of stored energy anymore (Ray, 1979). 
Furthermore, it remains unclear how to apply this method in case of 
complex circuits consisting of many loops. 



4-2 Problem 2: Self-Adjointness 

The reason for this discrepancy is that the differential 
equation (17) is not self-adjoint. To see this, let us consider 
the general form of a circuit consisting of memristors, me- 
minductors, memcapacitors, and their, possibly nonlinear, 
conventional counterparts, independent voltage sources, 
and independent current sources, given by 

Aij(x, x)x j + Bi(x, x) = 0, (19) 

where, in comparison to (8), it is observed that Aij now 
also depends on the x-coordinates. The necessary and 
sufficient conditions for the existence of a Lagrangian for 
this case extend to 

(20) 
(21) 
(22) 
(23) 



A, 



A, 
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> dAi ij.k 



dx k 



Returning to the differential equation (17), it is directly 
verified that, with A(q) = Lm{q) and B(q, q) = L' M (q)q 2 + 
q/C, the circuit is not self-adjoint, and therefore does not 
allow a Lagrangian formulation. 

A possible solution to compensate for the erroneous term 
is to add contra terms to the right-hand side of the 
differential equation; see Pesce (2003). This is tantamount 
to introducing a content function of the form 

V{q,q) = ^L' M {q)q\ 

like in (16). However, as argued before, the variational 
character is lost. 

4-3 Conservation of Flux and Charge 

The conventional Lagrangian formalism applied to electri- 
cal circuits essentially codes the KVL and KCL in terms 
of energy storage via the Lagrangian. On the other hand, 
suppose that we integrate both Kirchhoff laws with respect 
to time. This would result in the laws of conservation of 
flux and charge around a particular loop and node 

£>i(i)=0, f Vi(T)dT, (24) 

5>(*)=0, qj(t)= f 7,(r)dT, (25) 

respectively, which state that flux and charge can neither 
be created nor be destroyed (Chua, 2009). 

In the sequel, we will refer to (24) as the iKVL and to (25) 
as the iKCL. 

4-4 Memory State Functions 

In Section 2, we have seen that the most fundamental 
relationship of a charge-modulated meminductor is given 
by P = p(q), so that (6) is rather a consequence of the 



latter. It seems therefore more natural to consider, instead 
of the stored magnetic co-energy, which was obtained in 
(18) as the integral of (6) with respect to the current, a 
function of the form 



f*(q) 



p(q)dq. 



(26) 



When plotted in the g-versus-p plane, (26) represents 
the area above the curve associated to the constitutive 
relation. On the other hand, its complementary part, the 
function T, is defined to be the area below the curve and 
takes the form 

(27) 



f(p) := J q(p)dp. 



Under the assumption that p = p(q) and q = q(p) are 
invertible (i.e., one-to-one), the two state functions can be 
related via the Legendre transform 



T + T* 



QP, 



In a similar fashion, for a flux-modulated memcapacitor, 
we propose a function of the form 

(28) 



U*(ip) := J o-(tp)dip, 

representing the area above the curve associated to the 
constitutive relation, whereas for an integrated charge- 
modulated memcapacitor 

U(a) := J <p(a)da, (29) 

representing the area underneath the curve. 

Furthermore, we have 

Note that T, T* , U, and U* all serve as a state function and 
all exhibit the units of energy times time-squared (Joule- 
second-squared [Js 2 ]), which is equivalent to action times 
time. 

4-5 Meminductor and Memcapacitor Circuits 

Consider again the circuit of Fig. 2. First, we select, instead 
of a loop charge q, the integrated loop charge a, with a = 
q, as the configuration variable. Secondly, it is observed 
that the conventional capacitor can be considered as a 
linear memcapacitor with constitutive relation tp = a/C. 
In terms of the state functions proposed in the previous 
subsection, let us define a Lagrangian of the form 



f a 1 

C(a, a) = p{q)dq - —a 2 



Then, it can be demonstrated that, invoking Hamilton's 
principle of least action and considering variations in terms 
of (T, we get the Lagrangian type of equation 



d_0C _dL_ { , 
dt da da 



(30) 



which, in turn, generates the nonlinear differential equa- 
tion p'(a)a + a/C = 0. Differentiating the latter with 
respect to time yields 

p'{a)a+p"(a)a 2 +a/C = 0, 

which, after a change of variables a — q, and recognizing 
that the incremental meminductance L M (q) := p'(a), 



clearly coincides with the correct equation of motion (17). 
Furthermore, it is easily seen that (30) constitutes the 
iKVL for the circuit. Thus, we have rendered the dynamics 
self-adjoint by considering the circuit characteristics from 
the perspective of flux conservation. 

In general, the Lagrangian equations for circuits contain- 
ing charge-modulated meminductors, integrated charge- 
modulated memcapacitors, and their linear conventional 
counterparts, take the form 

dt d&i da- ~ U ' y6 ' 

with as Lagrangian £(cr, a) = T*(a) - U(a), where T*(&) 
and U(a) are now representing the sums of the individual 
memory storage functions associated to the meminductors 
and the memcapacitors in the circuit, respectively. 

Naturally, the dual version of (31) reads 

dt dpi dpi ~ ' 1 ' 

with the co-Lagrangian C*(p, p) = U*(p) - T(p). 

4-6 Including Memristors 

Although memristors are dynamical elements they also 
behave as resistors. For that reason they cannot included 
using a Lagrangian or co-Lagrangian alone, and therefore 
we need to introduce, at the cost of a underlying varia- 
tional principle, a second pair of state functions that play 
a similar role as the content and co-content functions for 
conventional resistors (see Subsection 3.2). 

For a charge-modulated memristor, we propose the state 
function 



<fi(q)dq, 



(33) 



which, when plotted in the y-versus-g plane, represents 
the area underneath the constitutive relation. 

For a flux-modulated memristor, we introduce the comple- 
mentary function 



q((p)dip, 



(34) 



representing the area above the curve. Note that (33) and 
(34) exhibit the units of energy times time (Joule-second 
[Js]), which is equivalent to action. 

As an illustration, consider the circuit depicted in Fig. 3 
consisting of a charge-modulated memristor Rm, a conven- 
tional linear inductor L, an integrated charge-modulated 
memcapacitor Cm, and a conventional linear resistor R. 
From a Lagrangian perspective, it is necessary to consider 
a loop analysis. For that, we select as configuration vari- 
ables the integrated loop charges cr l , with i = 1,2. The 
Lagrangian for the circuit reads 

1 _ 2 

C(a\a 2 ,a 1 ) = ~L(& 1 ) 2 -jT ' ^»da, 
whereas the memristive and resistive action is given by 

V(a\a 2 )=£ ^ M (g)d<z + ii?(a 2 ) 2 . 




0, i = 1,2, 



R 



Fig. 3. Circuit with a memristor and a memcapacitor. 

The associated Lagrangian equations take the form 

d dC dC dV 

dt d& 1 da 1 da 1 
which provide the iKVL for loop 1: 

La 1 + (pc M {a 1 - cj 2 ) + <p RM (a 1 ) = 0, 
and the iKVL for loop 2: 

-yc M {a l -cr 2 ) + Ra 2 = 0. 
The KVL's are obtained by differentiating the latter equa- 
tions with respect to time. 

5. ABSEMENT: A MECHANICAL ANALOGY OF 
INTEGRATED CHARGE 

In science and engineering, the ideas and concepts devel- 
oped in one branch of science and engineering are often 
transferred to other branches. One approach to transfer- 
ring these ideas and concepts is by the use of analogies. 
Classically, voltage is commonly considered as the electri- 
cal analogue of force, and current is the electrical analogue 
of velocity. Consequently, flux can be considered as the 
analogue of momentum or impulse, and charge as the 
analogue of position or displacement. 

In mechanics, displacement and its various derivatives 
define an ordered hierarchy of meaningful concepts. The 
first derivative of displacement is velocity, the second 
derivative is acceleration, the third derivative is jerk, 
the fourth derivative is jounce, etc.. On the other hand, 
recently also the integral of displacement over time is 
introduced to model the essential behavior of flow-based 
musical instruments, such as the hydraulophone discussed 
in Mann et al. (2006). This quantity is called absement, 
a contraction of absence and displacement, which, in SI 
units, is measured in meter times seconds [ms]. One meter- 
second corresponds to being absent one meter from an 
origin or other reference point for a duration of one second. 
Thus, integrated charge can be considered as the electrical 
analogue of absement. 

5.1 Absement Related to Memcapacitors 

To gain some intuition of absement in the context of the 
memory elements, consider the commonly used analogy of 
an electrical capacitor and a bucket that can be filled with 
water. If current were the flow rate of water from a tap, 
then charge would be the amount of water in the bucket. 
As the flow rate is proportional to how far open the tap is, 
the rate of flow increases as the tap is opened up further. 
Now, if we consider the displacement of the tap handle 
from its rest position, then the amount of water in the 
bucket is approximately proportional to the time-integral 



of the handle's displacement, i.e., the handle's absement, 
which is a measure of how 'absent' (how far and for how 
long) the handle is from its closed position. Equivalently, 
since displacement is the time-derivative of absement, the 
position of the handle is the time-derivative of how much 
water is accumulated in the bucket. 

From an electrical perspective, the absement of a capacitor 
is a measure for the amount of charge that is needed in a 
particular time interval to charge the capacitor to a certain 
level. For a memcapacitor this process involves both a 
charging and a shaping process of the capacitor and its 
capacitance. The amount of absement may then be used 
as a measure of how much charge and time is needed to 
write the memory of a memcapacitor. 

5.2 Absement Related to Memristors and Meminductors 

A mechanical example that exhibits meminductive phe- 
nomena concerns the elementary problem of a heavy cable 
that is deployed from a reel. As argued in Jeltsema and 
Doria-Cerezo (2010), under some reasonable assumptions 
on the geometry, this system admits a nonlinear constitu- 
tive relationship between integrated angular momentum p 
(the mechanical analogue of integrated flux) and angular 
displacement (the mechanical analogue of charge) a of the 
form p = p{9), with p — p. The amount of cable deployed 
from the reel equals R9, where R is the diameter of the 
reel. The absement in this case could be a measure for the 
amount of cable that needs to be deployed from the system 
in a particular time interval to empty the reel. 

The simplest physical example of a mechanical memristor 
is a tapered dashpot. This type of dashpot is a mechan- 
ical resistor whose resistance (i.e., its friction coefficient) 
depends on the relative displacement of its terminals. Al- 
though a physical electrical passive two-terminal memris- 
tive device was constructed only recently by Strukov et al. 
(2008), a tapered dashpot was already brought forward 
in the early seventies; see Jeltsema and Scherpen (2009). 
As for the meminductor, the amount of absement in this 
case could be a measure for the piston being absent from 
a certain resistance value for a duration of one second. 

6. FINAL REMARKS 

In this paper we have presented a novel (co-)Lagrangian 
framework to include memristors, meminductors, and 
memcapacitors, together with the conventional linear 
counterparts. In case of a circuit that consist only of me- 
minductors, memcapacitors, and their conventional linear 
counterparts, the Lagrangian equations can be obtained 
from Hamilton's principle of least action. Memristive and 
linear resistive elements can be included via the intro- 
duction of an action function that plays a role similar to 
the Rayleigh (co-)dissipation function or (co-)content. We 
conclude the paper with the following remarks. 

• Independent voltage (resp. current) sources can be in- 
cluded in the framework by considering them either as 
time-varying capacitors (resp. inductors) via a func- 
tion function of the form (29) (resp. (27)), or as time- 
varying resistors via a memristive action (resp. co- 
action) function (33) (resp. (34)). Furthermore, in 



order the bring the sources to the iKVL or iKCL level, 
we need to consider their time- integrals. For example, 
suppose we add a voltage source to the first loop of 
the circuit of Fig. 3. Then we need to extract either 
a term — a\(p e (t) from the Lagrangian or add a term 
—&\(p e (t) to the memristive action function, with 




where e denotes the source voltage. 

• The amount of integrated charge (electrical abse- 
ment) can be used as a measure of how much charge 
and time is needed to write the memory of a mem- 
ristor, meminductor, or memcapacitor. Although not 
addressed in the present paper, a similar interpreta- 
tion should apply to integrated flux, which in the 
mechanical domain corresponds to integrated mo- 
mentum measured in Newton times second-squared 
[Ns 2 ] or kilogram- meter [kg-m]. 

• Although we have restricted our attention to circuits 
that can be modeled via either a loop or a node analy- 
sis, mixed formulations are of course also possible; see 
e.g., Jeltsema and Scherpen (2009) for a discussion in 
the context of conventional circuits. 
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